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A MODEL FOR JET-NOISE ANALYSIS USING PRESSURE-GRADIENT 
CORRELATIONS ON AN IMAGINARY CONE 

By Thomas D. Norum 
Langley Research Center 

SUMMARY 

To determine the energy flux and far -field intensity of a jet through measurements 
of pressure-gradient correlations on a conical surface surrounding the jet, the appro- 
priate Green function is derived. The feasibility of this representation is checked by 
using a point source as the noise generator. Excellent agreement is obtained between 
exact results for the intensity normal to the cone and results from the Green function 
technique, which requires only a small amount of computing time. The representation of 
the far -field intensity at low frequency in terms of apparent sources on the cone shows 
the apparent sources to be spread almost symmetrically about the actual location of the 
point source as projected on the conical surface. The interpretation of this representa- 
tion at higher frequencies is less clear because of the presence of apparent sinks in the. 
source distribution. 


INTRODUCTION 

A method for determining the flux of acoustic energy from a jet by measurement of 
pressure -gradient correlations on an imaginary surface close to the jet is described in 
references 1 and 2. In this method the velocity field surrounding the jet is measured to 
establish the region where the velocity fluctuations satisfy the homogeneous wave equation. 
Once this acoustic region is determined, the imaginary surface over which correlations 
are to be measured is chosen to lie completely within this region, as close to the jet 
boundary as possible. At this inner limit of the acoustic field the pressure is entirely 
acoustical, and the measured pressure -gradient correlations can be related to the near- 
field and far -field acoustic energy through the wave equation. Hence a good estimate of 
the regions where noise leaves the jet and the contributions to the far field from these 
regions can be obtained from acoustic measurements alone. 

A plane surface near the jet was used for measurements of a subsonic jet in refer- 
ence 2. However, since the inner limits of the acoustic field of a circular jet form a 
conical surface enclosing the jet (fig. 1), the natural choice of a surface to give the most 
complete description of the radiation is a cone. 



The purpose of this report is to develop the mathematical model needed to deter- 
mine the near and far acoustic fields of a jet from pres sure -gradient correlations on a 
conical surface enclosing the jet. Although the use of this model is aimed at jet-noise 
experiments, its development is applicable to any source distribution enclosed by an 
imaginary surface. Hence a general formulation for the acoustical characteristics of the 
source in terms of the near-field correlation measurements is presented along with the 
necessary analytical developments for use with a conical surface. The intensity flux 
through the surface determined from this formulation is compared with the exact results 
for a known point source. The apparent source distribution on the surface (introduced 
in ref. 2) due to the known source is also obtained. These results validate the mathe- 
matical model and yield information concerning the interpretation of the apparent source 
distribution. 


SYMBOLS 


B 


i/m 


g 

G 


,( 2 ) 


H 

I 


defined in equation (54) 

frequency parameter, k g r s 

Green function constants (eq. (37)) 

ambient speed of sound 

unit vector in ith direction 

time -dependent Green function 

Green function in frequency domain 

spherical Bessel function of third kind of order v 

defined in equation (47) 

acoustic intensity 

spherical Bessel function of first kind of order v 
nondimensionalized intensity component normal to cone (eq. (43)) 
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wave number, a>/c 


integers 

normal vector to surface 

direction normal to surface measured positive into region U 
pressure 

associated Legendre function of first kind of degree v and order m 

apparent source distribution for far-field intensity 

associated Legendre function of second kind of degree v and order' 

position vector 

radial distance from origin 

distance; also separation function in radial direction 

real part of 

surface 

time 

averaging time 
sourceless region 
velocity 

spherical coordinate, cos 9 

0 


separation function in w-direction 



y 


p 


a, (3 


6 


e 


m 


e 


e 


o 


p 


p 

T 

<t> 

<P 

$ 

OJ 


V 

V 2 


spherical Bessel function of second kind of order v 

degrees of Legendre function; also nondimensional distance 

Dirac delta function 

1 if m = 0, 2 if m > 0 

zenith angle in spherical coordinates 

supplement of cone half-angle, see figures 1, 3, and 4 

degree of Legendre function 

density 

strength of point source 

separation time 

spectral-density function 

azimuth angle in spherical coordinates 

separation function in ^-direction 

correlation function 

angular frequency 

gradient operator 

Laplacian operator 


Subscripts: 

n direction normal to surface 

s pertaining to point source 
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0, 1, 2 distinguish variables over which integration is performed 
Superscripts: 

' Green function source coordinate 

* complex conjugate 

The notation ( ) denotes mean value. Arrows over symbols denote vector 
quantities. 


REVIEW OF GENERAL THEORY 


Since an expression is desired for the acoustical radiation from a jet in terms of 
near-field measurements on an imaginary surface near the jet, it is instructive to present 
the general formulation for the acoustical field in a sourceless region of space. Specific 
cases of this formulation have been given for the regions outside a cylinder (ref. 1) and 
opposite a plane surface (ref. 2). 

Let U denote this sourceless region which has as its inner boundary the surface S 
separating it from the region that contains acoustic sources (fig. 2). At any point r 
of U the pressure satisfies the wave equation 



±jA 

c 2 at 2 / 


p(r ,t) = 0 


Upon introducing a Green function in U which satisfies 



6(t-t') 


(1) 


( 2 ) 


and which satisfies an appropriate boundary condition on S, it can be shown by the pro- 
cedure of reference 3 that 


. \ r°° r [V \ ag(r ,F 0 ,t-t’) _ ap(rV T ) 

( r ’*) = L ^ Is dS ° P M 3n n ' ■ K" ’ r 0’ t ' t ’) “W 1 


( 3 ) 


where n denotes the direction normal to the surface S and is measured positive into 
the region U. 
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The boundary condition on g depends on whether measurements of pressure or 
pressure gradient are to be made. When the latter is chosen, the appropriate condition 
on g is 


8g(r,r ’,t-t’) _ 
an 


0 


(4) 


for all field points r located on the surface S. 

With t = t - t', the pressure field is then 

p ( F > 1 ) ' ' £«, dT £ ds o B ( F > f 'o- T ) ~ 


(5) 


The notation for correlation and spectral density is now introduced. For covariant 
stationary processes, the cross correlation of quantity A measured at the point 
with quantity B at point rL is denoted by 


T 

B ( r b- t+T )) Vfe^£ T dtA ‘( r a'‘) B ( f b" l+T ) 


( 6 ) 


The corresponding cross spectral density is 


^A,B( r a’V w ) = dT +A,b{Vv t ) e " ia,T (7 > 

If A is a vector quantity, the notation ij/~ g^a’^b’ T ) * s used ’ where 


( 8 ) 


i=x,y,z 


Subscripts x, y, and z denote directions along the Cartesian axes. The notation 
<p J B(^a , ^b’ w ) * s ex P resse ^ as in equation (7) in terms of g^a’^b’ 7 )' 

The correlation of pressure with velocity at the point F is 

?V,p( r ’ F ’ T ) = ' L d ‘ L 6S ° B ( r - F 0’‘) h 3£ ( F - F 0> T - 1 ) 


( 9 ) 
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Taking the Fourier transform of this result and manipulating the resulting integrals 
yield 


5v,pt r ' r '“) = - £ ds o 9p( r ’ r o>“) G ( r ’ f o-") 

V, 3n 

where G is the Fourier transform of the Green function g, that is, 

. r oo 

G (r,r 0 ,w).= j ^dre" la,T g(r,r 0 ,Tj 
From the equation of motion for the acoustic field, 



( 10 ) 


( 11 ) 


(12) 


the following relation between cross spectral densities is obtained: 


a( r - r o.") ‘ikKp SR ( r "’ r_ 0’“) 

. ’3n P, 9n 

Hence 

^ £ ds o ?_ a#-^) °( F - F o>“) 

F, 9n 


(13) 


(14) 


This relationship gives the spectral distribution of acoustic intensity at the point F 
in terms of the spectrum of the cross correlation of the pressure gradient at r with the 
pressure-gradient component normal to the surface S. Upon choosing r to be a point 
on S and taking the component of equation (14) normal to S, the spectrum of the normal 
acoustic intensity at F becomes 

■ *V n ,p( r < r >") - £ dS 0 % 3p( r > r 0- 1 ’) G ( F > F o>") <15 > 

3n’9n 

The spectrum is thus obtained in terms of the spectral distribution of, pressure-gradient 
correlations along the surface S. 
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The acoustic intensity normal to S is then obtained as 



(16) 


The distribution of radiated acoustic energy may be determined from the autocorre- 
lation of the far -field pressure. For any far -field point r in U, use of equation (5) 
gives this correlation as 

*P,p( F ’ F ’ T > ■ is dS l I s dS 2 II dt 2 II dt l 

xg*(r,r 1 ,t 1 ) g(r,r 2 ,t 2 ) ^ (F^F^T+tj-tg) (17) 

3n’an 

where r^ and r^ are points on S. Taking the Fourier transform of this relation and 
again manipulating the resulting integrals give 


Vp( F ’ r ■“} = £ dS l £ dS 2 G *( F °( F > T 2> a ) 4>ap 8 e( F X' F 2’“) (18) 


an’ an 


Hence the spectral distribution of the far -field intensity is also expressed in terms of the 
pressure -gradient correlations on S. The total radiated acoustic intensity at the 
point r is thus 


‘( r ) = £ ^p,p( r ' r ' °) “ 2 h~c L 0 p,p( ? ’ r ' 


(19) 


The far-field intensity is now represented in terms of an apparent source distribu- 
tion on the surface. This distribution determines the contribution to the far field due to 
apparent sources located on the incremental area dS of the surface and is defined by 


J ( F ) = I s ds i q i( r i ;F ) ■ (20) 

Combining equations (18) to (20) gives this distribution as 
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( 21 ) 


^(v) = 2^ L d “ tf ( F ’ F l’ a ’) Is dS 2 G ( F ' F 2.“) * 9p 2p( F l. r 2. w ) 

dn’ 9n • 

GREEN FUNCTION FOR A CONE 

To take advantage of the formulation given in the previous section, one needs to find 
the appropriate Green function. With t - t* = t equation (2) can be rewritten as 

g(r ,F, t ) = -6(r-r T ) 6 (t) (22) 

Taking the Fourier transform with k = u>/c yields 

(v 2 +k 2 ) G(r,r>) = -6(r,r’) (23) 

where for all points T on S, G must satisfy 

_8G(r ; r >) = Q (24) 

9n 

and must satisfy the condition for outgoing waves 

lim r { — + ikG) = 0 . (25) 

r—eo \ar . / • . ; .... 

The surface S is now chosen to be a cone enclosing the jet, as shown in figure 1. 

In terms of spherical coordinates (r,w ,<£), where w = cos 0 in the usual notation (fig. 3), 
equation (23) becomes 

(r-r 1 ) 6(w-w’) 6 (cp-cp') 
(26) 

The Green function G is constructed by considering the homogeneous form of this differ- 
ential equation. By the usual separation -of -variable technique, solutions of the form 

y(r,w,<p) = R(r) W(w) $(<p) (27) 
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are assumed, where y is a solution to the homogeneous form of equation (26). With the 
resulting separation constants set equal to m^ and v(v + 1), the following set of ordi- 
nary differential equations is obtained: 



(28) 


The three sets of linearly independent solutions to each of these equations are 


$((p) : cos(m<p), sin(m<p) 

W(w) ; P™(w), Q™(w) > 

R(r) : jj,(kr), y^(kr) > 


(29) 


where P™ and Q™ are the associated Legendre functions of the first and second kinds, 
respectively, and and y^ are the spherical Bessel functions of the first and second 
kinds, respectively (ref. 4). 

The Green function is now constructed from these solutions. Applying the condition 
of periodicity of cp over 2n shows that m must be an integer. Also, G must be 
symmetric with respect to therefore, the <p dependence is 


$(cp) oc cos 



(30) 


Since G must be finite at w = 1 (6 = 0), Q™ is eliminated from the solution 
because this function becomes infinite. The normal derivative on the cone is zero only if 
the degree v satisfies the condition 


dw 


w=w 


0 


= 0 


(31) 
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Hence 


W(w) cc P™( w ) (32) 

where p satisfies equation (31) and m is an integer. 

The condition of finiteness at r = 0 eliminates y^kr) from the solution for 
r < r'. At large values of r the outgoing wave condition requires 


R(r) a h^(kr) = j„(kr) - iy^kr) 

The solution that is continuous a cross the source point is thus 


(33) 


(kr) h( 2 )(kr') 




R(r) cc / 


j „(kr T ) h^(kr) 


(r £ r T ) 


(r s r') 






(34) 


J 


Combining the relations (30), (32), and (34) over all possible indices v and m 
gives the solution for the Green function as 

oo 

G = Y, I B um cos P ? (w) i^ kr< ) h l 2) ( kr> ) (35) 

m=0 v 

where m must be integral, v satisfies the condition specified in equation (31), and the 
notation 




F(x<) = 


F(x>) = 


F(x) 

(x £ x’) 

F(x’) 

(x ^ x') 

F(x') • 

(x 2 X’) 

F(x) 

(x S X') 


(36) 


J 


is used. 
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The constants B are obtained in the usual manner by substituting equation (35) 
into the differential equation (26) and successively: (a) multiplying by [cos dtp 

and integrating over ( 0,2n ); (b) multiplying by P^(w)dw and integrating over [wQ,lj , 

(the orthogonality of the associated Legendre functions being utilized as shown in the 
appendix); and (c) multiplying by r2 dr and integrating over the source point, that is, 
over.Xr'-Ojr’+O). Performing these calculations results in 

(37) 



Macdonald (ref. 5) shows that the values of v satisfying equation (31) form a dis- 
crete set of real numbers. Hence, upon utilizing the result (eq. (A7)) of the appendix, the 
Green function becomes - 


G = 


27ri ( 1 - w 0 2 ) nfe) 
(2v + 1) 


J e m cos 


P™(w’) P™(w) 


a dP™(wp 

P (w) ^ 

_ v dv dw 


jj,(kr<) h^ 2) (kr>) 


w=w„ 


(38) 


where the notation of equation (36) is used. Different representations for this function 
have been obtained by Carslaw (ref. 6) and Felsen (ref. 7) using more involved methods. 

CALCULATION OF GREEN FUNCTION 

The feasibility of using the conical Green function depends on fast and accurate cal- 
culation of Legendre functions of real degree and spherical Bessel functions of real order. 
The Legendre functions were determined via the Mehler-Dirichlet integral representation 
as described by Waterman (ref. 8). Results accurate to at least eight significant figures 
for the Legendre functions and their derivatives were obtained by this method. The 
degrees v needed in the Green function were determined so that the calculations in 
equation (38) are accurate to six significant digits. 

The Bessel functions were also calculated from an integral representation, with 
asymptotic forms used for both large argument and large order (ref. 4). Overlap between 
these methods of calculation assured an accuracy of six significant figures except near 
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the zeros of the Bessel functions, where accuracy to at least the sixth decimal place was 
obtained. - 

APPLICATION TO A POINT SOURCE 


To check the validity of the formulation and the practicability of its application with 
the minimum amount of mathematical detail, the jet was replaced by the simplest deter- 
ministic acoustic source, a point monopole of angular frequency w g located on the axis 
of the cone as shown in figure 4. The resulting intensity normal to the cone is computed 
by equations (15) and (16) and is compared with the result obtained by the direct method, 
which is described in the following paragraphs. 

The pressure field due to the point source is given by 




io> t -ik R 
i 5 e s 
47tR 


(39) 


-«* / — 2 — 2 

where R = r" - r^J = Jr g + 2wr g r + r and a g is the strength of the source. With 
the velocity field determined through the acoustic equation of motion (eq. (12)), the corre- 
lation of pressure and velocity becomes 


*p,v( r ’ r ’ T) = 


! ( k s R - 0 


VR e 


lw s T 


167T 2 pck R 3 


(40) 


Now choose r - to be located on the surface of the cone and note that 


9R 

an 




(41) 


along this surface; then the component of acoustic intensity normal to the cone is 


[ !i( r ) - R e (* Pl v n ( F ’ F ’°;j 


°s 2r s^ 


w , 


167T 2 pcR 3 


(42) 
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Nondimensionalizing this result as 



1677 3 pcr 3 

b 

"s 2 ^ 1 '* O 2 



(43) 


gives 



and hence the nondimensionalized intensity component normal to the cone is a function of 
only the relative distance along the surface rj r g) and of course the cone angle. 

To compare this result with that from the Green function approach, equations (15) 
and (16) are combined to give 


I„(?) - ^ 



(45) 


where r is any point on the cone. 

Calculating the pressure-gradient correlation by means of equation (39) yields 


^pap( r ’ r 0’ T ) 

0n’9n 


J s 2 ( 1 -”0 2 ) r s 2H ( k s B o) H *( k s R ) c to s T 
16tt 2 R 0 3 R 3 


(46) 


where 


H(x) s (1 + ix)e _ix 


(47) 


By using this expression the two infinite integrals are readily calculated; and nondimen- 
sionalizing again as was done for equation (43) yields 


J n(r) =Re- 


w 0 r s 


n* ( k s R ) £ H ( k s R o) „ 


R° 


dS 


Rn 




r (V W ( 


(48) 


If the Green function from equation (38) is introduced, and since an element of surface 
area of the cone is 
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where b g = k g r g . (The notation of eq. (36) is again used.) In contrast to the exact 
formulation (eq. (44)), J n here appears to depend also on the frequency through the 
parameter b g . 

The total radiated acoustic intensity at the far-field point r is given by equa- 
tion (19). For the point source being considered the pressure field from equation (39) 
yields 


I (r ) = s 

I dlrect ^ ' " le^pcr 2 


(51) 


since in the far field R = r. 

For this simple source it is of interest to determine the relationship expressed in 
equation (21), that is, the contribution to the far-field intensity from apparent sources 
located on the cone. Determining the intensity by means of equations (18) and (19) yields 


dS 2 G*( 


ap ap 1 

an’ an 


Substituting in the same way as was done to obtain the nondimensionalized intensity com- 
ponent normal to the cone (eq. (50)) gives 
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( 53 ) 


/ r \ 2 a 2 

’ (t) life 


V' * r°° r i H*(k Ft. ^ 

£v w) I 0 dri — ii yvj 


X 


Xv< 


dr 


r 2 H ( k s R 2), 


2 R^ 3 


V( k s r 2) 


where 


if— P°(w) 
A (w) = (2f + l)e 2 — - 

^ in 

9 dP “< W ) 


8v dw 


w=w. 


(54) 


and the asymptotic formula 


h t 2) ( k s r )' 


exp 


-UV - - 1 77 


k s r 


(as r — °°) 


has been used. In this expression the intensity appears to depend on the frequency of the 
source and on the far -field angle 6 = cos~^w, but from the direct expression (eq. (51)), 
it can be seen that this dependence is only apparent. 

Introducing for convenience the nondimensional coordinate /3 measured along the 
cone axis by 


ri 

0 = — cos = - 


w 0 r l 


(55) 


and normalizing the intensity by its value found by the direct method yield 


r- 1 — = y°° d 0 ^’(0; 
^direct ® 


w) 


(56) 


where 
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The real part of this expression is interpreted as the contribution to the far-field 
intensity at the angle 9 = cos“W from apparent sources located on the conical strip 
between /3 and /3 + d/3, and is analogous to equation (21) for the contribution to the far 
field due to apparent sources located on an incremental area of the surface. Its shape as 
a function of the two variables /3 and cu g gives an indication of how well the source can 
actually be represented by the apparent sources on the cone, and hence of how well these 
apparent sources can be retraced to determine an actual source location along the center 
line of the cone. 

Equation (57) also affords an additional check on the formulation, since an integration 
of the real part over all values of /3 should equal unity, whereas integration of the imag- 
inary part should yield zero. 

The limiting form of this function for low frequency is obtained by setting b g = 0. 
This yields 



RESULTS 

The point-source nondimensionalized intensity component normal to a cone of half- 
angle 12° is presented in figure 5. The solid line represents the direct values calculated 
from equation (44), whereas the plotted points were obtained from the Green function solu- 
tion (eq. (50)). These points were obtained for the frequency parameter b g equal to both 
1 and 10, (For jet-noise experiments, b g = 10 is estimated to correspond to a source 
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at the nozzle exit of frequency between 400 Hz and 4000 Hz, depending on the location of 
the cone vertex.) The number of eigenvalues needed to obtain a given degree of accuracy 
increased with frequency: Accuracy of 0.1 percent was obtained by using about 10 eigen- 
values for b s = 1 and 30 eigenvalues for b g = 10. The calculated points averaged 
about 1 minute of computing time on a CDC 6600 computer. 

The contributions to the far-field intensity at an angle of 45° to the cone axis due to 
apparent sources on the cone (eq. (57)) are shown in figure 6(a) for small values of the 
frequency parameter b s . The result for the case in which the frequency approaches 
zero (b g = 0, eq. (58)) is given along with results calculated for b g = 0.1 and 1.0. The 
distribution of apparent sources is spread almost symmetrically about the projected loca- 
tion of the actual source with little variation over this low frequency range. 

Figure 6(b) shows this distribution for higher frequencies of the source. The peak 
of the distribution again occurs at the projected position of the source, whereas the mag- 
nitude of the peak and the spread in the "downstream" direction increase with b s . 
Increasing frequency also results in negative contributions in the "upstream" direction, 
which correspond to apparent energy sinks in the interpretation given to equation (57). 

This result indicates that a high-frequency distributed source would appear to be located 
farther downstream than its actual location because of cancellation effects. Hence care 
must be exercised when making this apparent -source interpretation at the higher 
frequencies. 

The integration of the results for the real and imaginary parts of equation (57) does 
yield unity and zero, respectively, indicating accuracy of these results. The computing 
time for the far-field intensity was of the same order as for the near -field flux, the time 
increasing with the frequency of the source. 

CONCLUDING REMARKS 

The developments necessary for determination of the near and far acoustic fields 
of a jet in terms of measured pressure gradients on an imaginary conical surface sur- 
rounding the jet have been presented. A point monopole source was used in place of the 
jet to check the mathematical model. Accurate results for the intensity normal to the 
cone were found from this formulation in a small amount of computer time. At low source 
frequencies the representation of the far-field intensity in terms of apparent sources on 
the cone shows the sources to be spread almost symmetrically about the projected loca- 
tion of the actual source. At higher frequencies the apparent sources at the projected 
position of the actual source contribute an even greater amount to the far-field intensity, 
although upstream of the source the representation for the apparent sources becomes 
less clear because of the presence of apparent acoustical sinks. This phenomenon must 
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be recognized when the apparent source distribution calculated from pressure -gradient 
correlations is interpreted. 

The results indicate that the measurement of pressure-gradient correlations on a 
conical surface surrounding a jet is a practical method for determining the acoustic field 
of the jet. A more extensive analytic investigation should be made if an apparent-source 
interpretation of the far -field intensity is to be obtained from the experimental 
correlations. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 22, 1974. 



APPENDIX 


RELATIONS FOR ASSOCIATED LEGENDRE FUNCTIONS 


Let P™(w) and Pj^(w) be the associated Legendre functions of order m and 
arbitrary degrees a and /3. Hence they satisfy 




_d_ 

dw 


_d_ 

idw 


(l - w 2 ) ± 
dw 


(l - w 2 )— 
v 'dw 


Q!(Q!+1) - 


/3(/3+l) - 




1 - w 2 


m‘ 


1 - W 2 


>P“(w) = 0 


>P™(w) = 0 


(Al) 




Multiplying the first of these equations by P™(w) and the second by P^(w), subtracting 
and integrating w between Wq and 1 yield 

[pU 3 + 1) - a(a + 1)] J 1 dw P“(w) P™(w) 


0 


(l - w 2 ) 


dP^(w) 


dP^(w) 


P™(w) ^ - Pfi(w) ^ 

“ dw P dw 


1 w=w. 


(A2) 


Hence if the following condition is satisfied for v = a and /3: 


dP J “(w) 


dw 


= 0 


w=w. 


(A3) 


and if a ^ (3, then 


f 1 dw P>) P“(w) = 0 


(A4) 


that is, P™ and P™ are orthogonal on the interval ^Wq, lj . 

Since the degrees v of equation (A3) form a discrete set of real values (ref. 5), 
then for any root a there exists a value /3 = a + e, with e > 0, such that no roots exist 
in the interval (Qf,/3). Hence from equation (A2), 
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APPENDIX 


1 (l -W 2 ) P^(w)dP^ +£ (w) 
e 2a + 1 + e dw 


Therefore 



(A5) 


'-'W r 


dw 


P m (w) 
a ' ' 


r l 

= lim \ dw P m (w) P m (w) 

£-0J w . “ «+e 


(l -w2) P“(w) Um l dP “ +E ( w > 


2a + 1 


e—0 6 dw 


(A6) 


W=W r 


But the limit in this last expression is simply the rate of change of dP™(w)ydw with 
respect to the degree a. Hence 


I 1 

J w r 


dw 


P >> 


2 


(l - W 2) P™(w) 8 

dP“(w) 

2a + 1 da 

dw 


(A7) 


w=w r 


for any a satisfying equation (A3). 
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